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ABSTRACT 


Liquid crystals, a material that exhibits brilliant changes in 
color over narrow temperature bands, have been successfully used to 
study the temperature field produced by spherical and cylindrical 
radio frequency surgical probes. An agar-water mixture was used to 
simulate tissue. Experimental temperature data were obtained by 
photographing the liquid crystal material through the clear agar. 
Comparison of experimental data with an existing theoretical solution 
was excellent, agreeing within 8% for the spherical probe. However, 
the cylindrical data and cylindrical theory agreed only to within 30%, 


the large discrepancy being attributed to edge effects. 





~ TABLE OF CONTENTS 


i PRODUCTION --22225---_-e5----.-.--.-.-----2 
Il. HISHORMeAUMEAOKGROUNDM@ees........2.-----...__- 
Me. PHEOREPICAL DEVEWOPMENT® _....._...........____- 
py SPHERTCHEEPROBE MEE nema. eee ee 
B. CYLINDRICAL PROBE --.---------------------- 
IV. eS Ns es, rr 
io EQUBPNENT. ...csieee. Ste. oscees. eee. 
B. EXPERIMENTAL PROCEDURE .__._-_-_-_-.--_-__- 
V. BAPEREMENTALL RESULTS -~s8.------..ewe.-..-.-..- 
fer GE RIe MORO EES. Seekers... 205 22 
SE a Cc ae 
ml D.LSCUSGIGNMOF ARES SHeeeeteeteee._-_____seeeewe 
VII. CONCLUSIONS AND RECOMMENDATIONS -----..------.- 
i CONC SNS = i ne 
SP ECOL 0 Sr 
Pree Ac RESIS WIV UI Ngee one oe nee 
PePEWOIX B: LIQUID CRYSTALS ----ceo--.-.........-2..-.- 
eC: FUMED LORS 2a 8 ee 


PmeeNDIX Diss SAMPLE UNCERTAINTY ANALYSIS ..-_-.....---.-. 
PoeereiX Eb: COMPUTER PROGRAMS ..._....................- 


APPENDIX F: EXPERIMENTAL DATA 


a ae 


INITIAL DISTRIBUTION LIST 
FORM DD 1473 





i. 


2. 


is. 


LIST OF ILLUSTRATIONS 


Richardson's Technique applied to a spherical probe in a 


Per me ace cima le mnaS)St1V LY -----aea---—-- =~ ~~~ ----- === == 21 
Experimental equipment --------------------------------------- 24 
Circuitry for radio frequency probe experiments -------------- La 
Probe test celIS ----------n- nnn nen enn nnn nnn nen nnn nn nen ee 26 
Liquid crystal holders --------------------------------------- 28 


. Spherical transient temperature field growth with probe 


diameter of 0.241 cm, ------------------- SERS SSS SSS SS S18 
Nondimensional isotherm radii vs. nondimensional time for 

a spherical probe of diameter 0.482 cm. ---------------------- 34 
Nondimensional isotherm radii vs. nondimensional time for 

a spherical probe with diameter 0.482 cm. using three 

liquid crystals a a ia i pS 35 
Cylindrical transient temperature field growth 

L/D = 4, probe diameter = 0.241 cm. -------------------------- 3/ 
Cylindrical transient temperature field growth 

L/D = 10, probe diameter = 0.241 cm, ------------------------- 38 
Cylindrical transient temperature field growth 

L/D = 20, probe diameter = 0.241 cm, ------------------------- 39 
Nondimensional isotherm radii vs. nondimensional time 

for a cylindrical probe with L/D = 4 and diameter 0.241 cm.--- 40 
Nondimensional isotherm radii vs. nondimensional time 


for a cylindrical probe with L/D = 10 and diameter = 0.241 cm- 41 





14. 


nD. 


a 


HO . 


ie. 


ae. 


a. 


ae. 


28. 


24, 


ZS. 
20. 


Nondimensional isotherm radii vs. nondimensional time 

for a cylindrical probe with L/D = 20 and 

coerce — OFC 4M). .222--2-- one eeee--- ae a — 
Liquid crystals map a nonsymmetric thermal region ---------- 


Constant current lines emanating from a cylindrical 


Comparison of experimental data with original theory 

and corrected theory for cylindrical probe ----------------- 
Liquid crystals demonstrating the pronounced edge 

effects caused by a cylindrical probe ---------------------- 
Boiling in the agar-water mixture -------------------------- 
Color photograph of liquid crystals in thermal field 

of cylindrical probe (L/D = 4). Current 

a OOM meade CMe wee oO —— 22S -_. .-_------.---.--- 
Color photograph of liquid crystals in thermal field 

of cylindrical probe (L/D = 10). Gurrent 

flux = 0.0414 amps/cm? ----------------------------n en en eH 
Color photograph of liquid crystals in thermal field 

of cylindrical probe (L/D = 20). Current 

flux = 0.0414 amps/cm¢ -------------------------------------- 
Color photograph of liquid crystals in thermal field 

of spherical probe. ----------------------------------------- 
Comparison of cylindrical data and theory with spherical 
theory and scaled cylindrical data -------------------------- 
The resistivity test cell ----------------------------------- 


MMNALIVity teat call circuitry ----.--2---.--_----------_ == 





Gl. 


28. 


Zo. 


30. 


Sil. 


Bc. 


e3. 


34. 


5. 


56. 


Gi. 


So. 


ce 


ert by best Cell In the OVEN "TTT ttt Tt tt onto nnn 62 


Variation in resistivity values ----------------------------- 63 
Nondimensional resistivity vs. nondimensional temperature --- 65 
Birefringence and circular dichroism of liquid crystals ----- 68 
eee oe eee ge = 75 


Isopotentials and constant current lines for a 

cylindrical probe (L/D = 4) --------------------------------- 76 
Isopotentials and constant current lines for a 

cylindrical probe (L/D = 10) ------------------------------+- 17 
Isopotentials and constant current lines for a 

Sucndeaueprowe (Ween 20) 222-..--..-..-i---.ssaeeee=-.-.- 78 
Isopotentials and constant current lines for a 

cylindrical probe (L/D = 4) with an additional 

conductive boundary ----------------------------------------- 80 
Tsopotentials and constant current lines for a 

cylindrical probe (L/D = 10) with an additional 

Conguemi emeotndalit(e--———-.--—.—-—-----—-_ ~~ ------ - =~ - = - = - === 8] 
Isopotentials and constant current lines for a 

cylindrical probe (L/D = 20) with an additional 

conductive OE Bis eee a ee 82 
Isopotentials and constant current lines for a 

cylindrical probe in a plane perpendicular to the 

DYODE AXIS wen nen nen nr nnn rn ene en ne nnn nn nnn ener enc nn nae 83 
Tsopotentials and constant current lines for a 

cylindrical probe in a plane perpendicular to the 


ODO CG cela ae ne et ama a 83 





Symbol 


o om >o YP 


10 


NOMENCLATURE 


Description Typical Units 
Surface area cm2 
Wavelength Angstrom 


Number of flow tubes 
Volumetric heat capacity cal/cm? °C 


Temperature in degrees Centigrade 


Diameter cm 
Electric potential volts 
Frequency hertz 


Nondimensional power 

Corrected nondimensional power 

Current density amps /cm? 
Radio frequency current amps 
Coordinate of nodal point along 

distance axis in finite 

difference formulation 

Thermal conductivity W/m °K 
Length cm 


Stability criterion for finite 
difference formulation 


Slope of nondimensional resistivity 
curve 


Coordinate of nodal point along 

time axis in finite difference 
formulation 

Nondimensional radius to an isotherm 


Ten ohm-ten watt resistor 





Probe radius cm 


aS 
r Radius cm 
re Isotherm radius cm 
T Temperature ne 
i Time seconds 
V Voltage volts 
0 Thermal diffusivity cm?/sec 
Y Slope of dimensional resistivity 
curve 
n Ratio of speed of light in air 
to that in water 
6 Nondimensional temperature 
A Wavelength | Angstrom 
0 Resistivity | ohm-cm 
6 Nondimensional resistivity 
T Nondimensional time 
> General potential 
W Uncertainty 
V Del operator 
SUBSCRIPTS 
O Initial conditions 
10 At the ten ohm resistor 
Gat In the circuit 
7 Of the ith component 
D From the photographs 
cyl With reference to the cylindrical probe 
sphere With reference to the spherical probe 
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I. INTRODUCTION 


In neurosurgery, radio frequency current has been used to create 
lesions, thus alleviating conditions such as acute pain and Parkinsonism. 
A current emitting electrode, similar to a hypodermic needle, is inserted 
into the brain through the skull]. An inactive electrode is placed else- 
where on the body in an area of high conductivity. Radio frequency 
current, at a low power level, is passed between the electrodes and 
electrical energy is transformed to heat in the tissue with the result- 
ant rise in temperature destroying the tissue. 

To accurately use the radio qaeatietar technique, a surgeon must 
know the exact location of the electrode (probe) tip in the brain, the 
volume of tissue to be destroyed, the thermal and electrical properties 
of the tissue [1], and the temperature near the probe tip. Thermocouples 
and thermistors attached to the probe have been used to monitor probe tip 
temperatures [2, 3, 4]. However, these temperature sensors indicated the 
temperature of the probe, not the tissue temperature itself. Addition- 
ally, the sensors were nearly as large as the probe tip and may have 
distorted the thermal -field by acting as thermal sinks. Without accurate 
Knowledge of tissue temperature, tissue boiling occurred with attendant 
gas formation and tissue carbonization [4]. These effects resulted in 
irregular shaped lesions and unpredictable damage along the entrance 
path of the probe into the tissue. 

The objectives of this thesis were twofold. First, to experimentally 
determine the thermal field produced by a radio frequency probe embedded 
in a resistive agar medium. Liquid crystals, a material that indicates 


temperature through a color change, were used as the temperature sensors. 


1] 





The second objective was to compare experimental data with results pre- 
dicted using a theoretical heat transfer model of the radio frequency 
probe [5]. The general theoretical model included the effects of tissue 
metabolism and blood profusion. These two effects have been excluded 
from the present investigation since the thermal field was generated in 
an inert gel-like substance. 

The temperature fields produced by spherical and cylindrical radio 
frequency probes were successfully studied with the liquid crystal 
material. , The crystals, with their low thermal capacity [6], produced 
brilliant, undistorted, continuous displays of several selected isotherms 
Surrounding the probes. Agreement between the theoretically predicted 
and experimentally measured values of temperature for the spherical 
probe was excellent (8%). Agreement between theory and experiment for 
the cylindrical probe was not good (30%). The large differences observed 
in the cylindrical case are attributed to edge effects. The theory assumed 
an infinitely long probe while the experiment employed probes with length 
to diameter ratios of 4, 10, and 20. The edge effects were clearly dis- 


played by the liquid crystals. 
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II. HISTORICAL BACKGROUND 


Use of radio frequency current to create physiologic lesions is not 
new. Carpenter and Whittier [7] describe early investigations into the 
uses of radio frequency current as a surgical tool. Aronow [1] was the 
first investigator of recent time to study the problem of radio frequency 
current in tissue. Lesion formation using radio frequency current is 
accomplished by inserting into the tissue a small hypodermic needle which 
is electrically insulated except at the tip. This electrode acts as the 
active, current emitting electrode. A larger, inactive electrode is 
placed elsewhere on the body. The path of current conduction, i.e. the 
body, acts as a leaky dielectric by dissipating electrical energy in the 
form of heat [1]. This Joulean heating effect results in a local rise in 
tissue temperature. When tissue temperature is increased to 55°C, tissue 
destruction occurs. Accurate prediction of when tissue in the probe tip 
region has reached this temperature has been a major problem with the 
radio frequency probe. 

Excessive tip temperature causes gas formation and tissue carboniza- 
tion which results in poor prediction of lesion size and additional 
damage along the path of electrode insertion [4]. Studies with thermo- 
couples and thermistors have been conducted to correlate electrode 
temperature with lesion size [2, 3]. Since heating takes place in the 
tissue and not in the electrode, the sensors inside the probe tip indi- 
cate only that temperature rise caused by conduction from the tissue back 
into the electrode. This temperature rise may be less than the maximum 
tissue temperature. Therefore, reliable prevention of excessive tip 


temperature has been difficult. The same type of tissue temperature 


measurement problem occurs with the direct current probe, [5]. 
13 





The direct current probe has a more serious liability than the 
tissue temperature measurement problem. Direct current has preferred 
paths in the body which results in irregular lesions. Death may occur 
if the direct current finds a path to the heart. Preferred paths in 
the body do not appear to be a problem with the radio frequency or the 
_ high resistance probe. | 

The high resistance probe consists of a high resistance iron-nickel 
alloy wrapped around a nonconductive rod. Electrical current is passed 
through the wire causing the probe to be heated. Groff's study of the 
high resistance probe [8] demonstrated the feasibility of producing 
controlled, predictable heated regions in agar. In opposition to the 
radio frequency probe, where heating occurs in the tissue, the high 
resistance probe itself generates heat. Thermal energy is then trans- 
ferred to the tissue by conduction. Prediction of the magnitude of the 
temperature field produced by a high resistance probe thus requires 
knowledge of only the thermal properties of the medium. Another advan- 
tage of the high resistance probe is the fact that a surgeon may check 
for proper probe tip location before forming the lesion. 

Probe tip location may be verified by raising the tissue temperature 
to 40-49°C using low current levels. Increasing temperature into this 
range blocks nerve functioning but does not cause permanent damage. This 
technique, equally applicable to the high resistance and the radio fre- 
quency probe, is especially helpful in treating Parkinsonism. As the 
temperature increases, shaking diminishes and stops when the probe tip 
has been properly located. Lesion formation is then accomplished by 
increasing the local tissue temperature to 55°C, being careful to prevent 


excessive tip temperatures that result in tissue boiling. 
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While the high resistance probe requires knowledge of only the 
thermal properties of tissue, the radio frequency probe requires 
knowledge of both thermal and electrical properties. Geddes and 
Baker [9] discuss difficulties in obtaining electrical properties of 
tissue. Since tissues are composed of cells, different properties 
may be exhibited depending on the cell orientation. Other factors 
influencing the electrical property values include age, history, and 
environment of the tissue. Tissue properties measured IN VITRO 
differ from those measured IN VIVO. Thermal properties, on the other 
hand, are reliably predicted based on the tissue water content [10]. 
The prediction of lesion size based on the magnitude and duration of 
radio frequency current and tissue electrical properties is, there- 
fore, more complex when using the radio frequency probe than when 


using the high resistance probe. 





III. THEORETICAL DEVELOPMENT 


A. SPHERICAL PROBE 

Consider a spherical probe, of radius ie which continuously emits 
radio frequency current, I, into a medium that is homogeneous and infinite 
in extent. Assume that the medium possesses electrical resistivity, p, 
thermal conductivity, k, and volumetric heat capacity, C. Gengler [5] 
has derived the basic equation which accounts for the thermal effects 
produced in the medium due to the Joulean heating effects of the radio 


frequency current: 


r- ar 1672r att 


| 37 | 

ils 2. 

lca alee ¢ al (1) 
In the present investigation, the experimental test medium, an agar- 


water mixture, was found to exhtbit an electrical resistivity that 


varied linearly with temperature: 
Poets Alcan (2) 


where the subscript indicates an initial, reference condition. 


Substitution of Equation (2) into Equation (1) yields: 


Q 
— 


oT 
1 dace 3 , Izy (T - To) | 
z res i + = OF (3) 


Assuming that k and C are constant, Equation (3) may be normalized by 


introducing the following non-dimensional quantities: 
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R = a Distance 
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m - — 
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where I = coe = * the current flux. 


Substitution into Equation (3) yields: 


2 20] 
al RL , oy 4) = 38 (4) 
R2 oR Ro T 


The boundary and initial conditions are: 


a) at the probe surface, there is no heat conduction back into 


the probe, 
30 = 0; (5) 
dR Pe 


b) the temperature field is undistorted at large distances 


from the. probe, 


R+«,@e@+ 0; (6) 








c) the medium is initially at uniform temperature, 


—= Osea — 0. 


(7) 


Equations (4), (5), (6), and (7) may be cast in explicit finite dif- 


ference form as [5] 
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(iaR)* At 
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where i and n define the coordinates of a nodal point in an R, t 
network such that the nondimensional temperature at time t = NnAt 
and radial location R = iAR is defined as 04 s 17 1S ae KR = Ps 


i= i*, i*41,...,N-7,N. Solving Equation (4a) for of” 


yields: 


(4a) 


(6a) 


(7a) 
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Since an explicit finite difference scheme has been employed, a stability 
criterion must be stated to limit the size of the time step based on the 
value of incremental distance (AR). Discussions of stability may be found 


in References 11 and 12. The stability criterion for this problem becomes 


MmR2G 
O 


Solving for M and assuming the worse case, 


I (9) 


mAR2G 
TAR)* 


2 . 
Since M was originally defined as At/AR2, choosing a value of AR? fixes 
the value of At. 

Computer program I (Appendix E) was used to solve Equation (8) subject 
to boundary, initial, and stability conditions (Equations 5a, 6a, 7a, and 
9). Reference 13 states that the error in an explicit difference approxi- 
mation is on the order of AR2. By varying the size of the AR increments 


and plotting temperature versus AR?, it was possible to extrapolate the 


exact value of temperature. Figure 1 represents this technique, known as 
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Richardson's technique, for the spherical probe. It can be seen 
that using a value of AR = 0.1 yielded theoretical results within 


2% of the exact value. 


B. CYLINDRICAL PROBE 
Paralleling the development of the basic equations for the spheri- 


cal probe, the equation for the cylindrical probe in nondimensional 


fromm 1s: 
06 
i 0 R aR A G. i 38 110) 
R aR R2 It 
2 2 
where G = PaO ats 
kT 
O 
m= 1 = I 
A eur L 
O 


Substituting G = 6 + me) into Equation (10) and transforming to 


finite difference form yields: 


GM GM 
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The appropriate boundary conditions are: 
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where 1* corresponds to R = 1.0 
Moe We ow Is..., N=1, XN 


The stability criterion for the cylindrical case is: 


where M is defined as in the spherical case as At/AR?. 
Computer program II (Appendix E) was used to generate the numerical 


solution to this problem. 
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IV. EXPERIMENT 


A. EQUIPMENT 

A radio frequency generator, Model RFG-2AV, manufactured by 
Radionics Incorporated of Burlington, Massachusetts, was used in the 
present investigation to generate radio frequency current. Output 
from the generator was a continuous, unmodulated sine wave at 500 
kilohertz with a maximum power of 18 watts. Although the generator 
included a built-in radio frequency voltmeter and radio frequency 
milliammeter, two independent RMS meters were utilized in data collec- 
tion to allow the use of more sensitive scales where required. Figure 2 
displays the generator, voltmeters, and one of the experimental probe 
test cells used in the experiments. Frgqure’3 is a schematic of the 
circuitry used. 

Both spherical and cylindrical radio frequency probe geometries 
were Studied experimentally. An agar-water mixture (a gel-like 
substance) was used to simulate tissue. The mixture was contained in 
the test cells shown in Figure 4. In both cells, at least 25 probe 
radii separated the probe from the nearest solid boundary. Thus the 
cells could be considered as surfaces at infinity in relation to the 
probes. Plexiglass was used for the probe-agar cells due to its trans- 
parence and because it is a dielectric at the low power levels used 
during experimentation. Boundary conditions were changed from insulated 
surfaces to conductive surfaces by lining the inside surfaces of the 
cells with conductive foil. The current path was then from the probe, 
through the agar medium, to the foil with electrical energy dissipated 


as heat in the aqar. 
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PROBE CELL RADIO FREQUENCY 
GENERATOR 


Voltmeter Vy measured the voltage drop across the 
ten ohm resistor Ra9 . Circuit current was deter- 
mined from the two quantities using Kirchoff's Law. 


Cireuaury for radio frequency probe experiments 


Figure 3 
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Figure 4. 








Local temperature increases in the agar were measured by use of 
liquid crystals. Appendix B presents a complete discussion of liquid 
crystals. Briefly, liquid crystals change color at specific temperatures, 
the particular temperature depending on the crystal in use. The color 
change thus allows a continuous display of isotherm fronts for the map- 
ping of thermal regions. 

In the experiments, liquid crystals manufactured by National Cash 
Register and Hoffman-LaRoche were used. The Hoffman-LaRoche liquid 
crystal, ROCHROME, was provided in the form of a tape which could be 
applied directly to the object under study. The National Cash Register 
liquid crystals were obtained in. liquid form. The crystals are actually 
encapsulated in spheres 5 to 30 microns in diameter, thus preventing 
contamination of the crystal. Due to this encapsulation procedure, it 
was possible to mix different liquid crystals together without changing 
the characteristics of any single crystal. Several distinct sets of 
isotherms were viewed since the liquid crystals mixed did not have tem- 
perature ranges that overlapped. 

In the experiments, three liquid crystals were mixed together and 
were applied to a transparent 0.003 inch mylar substrate. The liquid 
crystals were then covered with black paint and polyurathane varnish. 
The black background was required to absorb the light transmitted 
through the liquid crystals; the polyurathane varnish sealed the black 
paint. Sealing the system in this manner allowed repeated use of the 
liquid crystal system. Figure 5 displays the liquid crystal sheet 
encased in the holders used in the test cells. 

Prior to experimental trials, the crystals were calibrated using 


the Rosemount Constant Temperature Bath following the procedure of 


er 
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Petrovic [14]. The following liquid crystals were calibrated and 


subsequently used in experimental trials: 


Table 1 
Manufacturer Designation Red Green Blue 
National Cash Register R-2/ C67 G ZO G 28.9°C 
Hoffman-LaRoche Rochrome 259) C Sino oe 34.2°C 
National Cash Register R-3/ Shae S720°€ ete peor & 
National Cash Register R-53 S065 Bec Saeo° 


Data were collected at specified times after application of power by 
photographing the isotherms displayed on the liquid crystal sheet. 
Distances from the probe center to an isotherm (color band) under consid- 


eration were easily measured from the photographs. 


B. EXPERIMENTAL PROCEDURE 

The agar solution was prepared by adding 3.0 grams of DIFCO 
“BACTO-AGAR" powder to 997 grams of boiling distilled water. When the 
powdered agar had thoroughly Beeoived: a 100 ml. sample was taken to 
fill the resistivity test cell which was used to determine the variation 
of resistivity with temperature (see Appendix A). Prior to jelling, the 
remaining agar solution was poured into one of the probe test cells. The 
liquid crystal sheet was positioned and the probe located such that a probe 
diameter was in the plane of the liquid crystals. The entire apparatus was 
then cooled for approximately six hours to achieve a uniform temperature in 
the agar. At that time the initial agar temperature was recorded, the 
generator circuit was connected, and power was applied. Photographs of 
the liquid crystal sheet were taken through the agar at specified times, 


thus recording the color change in the liquid crystals. Distances were 
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measured directly from the pictures and were scaled with respect to 


the probe radius to form a nondimensional distance. 


All quantities were nondimensionalized to allow general application 


of the data. 


De 


where 


The following nondimensional quantities were used (5): 





i 
- Distance 
O 
= Time 
O 
a 
7 S Temperature 





2 
Ite 6 Power 
A KT | 


1 + me Resistivity 


surface area of the probe (cm?) 

for the sphere Any * 

for the cylinder én ob 

distance from the probe center to the outer edge of the 
isotherm (color band) under consideration (cm) 
probe radius (cm) 

thermal diffusivity (cm2/sec) 

time (sec) 

isotherm temperature (°C) 

initial agar temperature (°C) 

radio frequency current (amps) 

initial agar resistivity (ohm-cm) 


thermal conductivity (W/cm °K) 
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Using a constant temperature-power ratio (9/@ ), curyes were 
plotted for distance (R) versus time (+) based on the solution of 
the finite difference equations for the probe geometry under investi- 
gation (Equations 8 or 11). The curves were then compared to the 
experimental data points of nondimensional isotherm location (R) at 
various times (t). An uncertainty analysis was conducted and is 


presented in Appendix D. 


Sil 





Y. EXPERIMENTAL RESULTS 


A. SPHERICAL PROBE 

Figure 6 is a series of photographs of the growth of the radial 
temperature field around the spherical probe at nondimensional times 
of 1.5, 3.1, 6.2, and 12.5 (corresponding to times of 60, 120, 240, 
and 480 seconds for a spherical probe of radius 0.24 cm.). In 
Figure 6c, the isotherms used in data collection are labeled for 
reference. The temperature-color relationships used in data 


collection were: 


Liquid Crystal Temperature Color 
R-27 | 7 oat Green 
Rochrome 2-9 20 Red 
R-37 So7ae C Green 


The R-53 appeared on the photograph at large times, but was not 

actually used in data reduction and comparison with theory. Values of 
nondimensional radius (R) and time (1+) were plotted for comparison with 
the theoretical curves based on temperature-power ratios (0/G ). 

Figures 7 and 8 are representative of the comparisons. The dotted 

lines are the uncertainty limits, which are discussed in Appendix D. 
Experimental results were within 8% of the theoretical values. Thus, the 
Spherical theory, taking into account variable resistivity, was verified 
by the experimental procedure of using liquid crystals to determine 


thermal changes. 
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Figure 6a Figure 6b 
Tau = 1.5 t = 60 sec. Tau = 3.1 t 38120 sec, 


Green BeZ7| 


Red Rochrome 





Figure 6c Fig 
Tau®="6.2 t = 2408sec. Tau = 12.5 


e 6d 
t = 480 sec. 


Figure 6. Spherical transient temperature field growth 
with probe diameter of 0.241 cm. 


33 





Figure 7 


Nondimensional isotherm radii 
VS. 
nondimensional time for a spherical 
probe of diameter 0.482 cm. 
Data set I 





4.0 
5 
"i 
a 
oO 
2.0 
Theoretical] curve 
i. Uncertainty limits 
© Experimental data 
0 
0 5 10 15 20 
Tau = a 





Figure 8 


Nondimensional isotherm radii 
VS. 
nondimensional time for a spherical 
probe of diameter 0.482 cm 
using three liquid crystals. 
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B. CYLINDRICAL PROBE | 
Figures 9, 10, and 11 display the experimentally determined 

transient temperature field at nondimensional times of 4.2, 6.25, 
9.3, and 12.5 (corresponding to times of 40, 60, 90, and 120 seconds 
for a cylindrical probe of radius 0.12 cm.). The length to diameter 
ratios used were 4, 10, and 20 in Figures 9, 10, and 11, respectively. 
Values of distance (R) were plotted as a function of time (t) for 
constant temperature-power ratios (0/G_), as seen in Figures 12, 13, 


and 14. The solid lines represent the theoretical curves. 
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Figure 9a Figure 9b 
fa = 4.2 t = 40 sec. fau= 6.2 t =%0'sec. 





Figure 9c Figure 9d 
Tau = 9.3 t = 90 sec. Tau = 12.5 t = 120 sec. 


Figure 9. Cylindrical transient temperature field 
growth, L/D = 4, probe diameter = 0.241 cm. 
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Figure 10a Figure 10b 
Tau = 4.2 #42t = 40 sec. Tau = 6.2 +t = 60 sec. 





Figure 10c Figure 10d 
Tau = 9.3 t = 90 Sec. Tau = 12.5 t = 120 sec. 


Figure 10. Cylindrical transient temperature field 
growth, L/D = 10, probe diameter = 0.24] cm. 
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g 
52 t= 607sec. 





g Figure lid 
3 t = 90 sec. Tau = 12.5 t = 120 sec. 





Tau = 


Figure 11. Cylindrical transient temperature field 
growth, L/D = 20, probe diameter = 0.24] cm. 
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Figure 12 
Nondimensional isotherm radii 


VS. 
nondimensional time for a cylindrical 
probe with L/D = 4 and diameter = 0.241 cm. 
Data set M 
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Figure 13 
; Nondimensional isotherm radii 
10.0 a 
nondimensional time for a cylindrical 
probe with L/D = 10 and diameter = 0.241.cm. 
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Figure 14 
Nondimensional isotherm radii 
| VS. 
nondimensional time for a cylindrical 
probe with L/D = 20 and diameter = 0.241 cm. 
Data set K 
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YI. DISCUSSION OF RESULTS 


A serious difficulty in modeling the radio frequency probe in 
tissue has previously been an inability to accurately map the thermal 
area surrounding the probe. Thermocouples and thermistors have either 
been placed inside the probe tip, where they actually measured the 
probe temperature, or placed in the medium, where they acted as thermal 
Sinks. The use of liquid crystals, as demonstrated in this thesis, 
now allows an accurate, continuous display of specific temperatures 
depending on the crystal in use. Liquid crystals provide a flexible 
tool for temperature measurement regardless of probe geometry. Figure 15 
Shows liquid crystals indicating isotherms which were caused by radio 
frequency current between two nonsymmetric probes. As can be seen, the 
liquid crystals react to provide a complete map of the thermal region. 

In the spherical geometry, the experimental values were in good 
agreement with the theory as modified for variable resistivity (8%). 
The cylindrical data and theory were not in agreement (30%). As a result 
of the poor agreement between heer and experiment in the cylindrical 
geometry, an examination of the flux lines emanating from the probe 
was conducted. Appendix C discusses the flux plots obtained. 

Figure 16 is a flux plot showing only the lines of constant current. 
As can be seen, the direction of the current lines is changed in the 
area above the probe. The bending of the current lines results ina 
higher local current density with greater heating than theoretically 


predicted. 
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Figure 15 
Liquid crystals map a 
nonsymmetric thermal region 
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Insulated Boundary 
(Top of Agar-water mixture) 


es 
6. 
: 
Line of Conductive 
Symmetry Boundary 





Conductive Boundary 


Figure 16. Constant current lines emanating from 
a cylindrical probe based on teledeltus plot of 
the isopotentials in plane perpendicular to probe 
axis 
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An approximation was derived to account for the increased current 
density. Due to symmetry of the problem, only half the probe-cel] 
system was modeled. The following parameters were defined: 

Number of flow tubes (for half field) = b 

Current flow (for half field) = 1/2 

Half the probe surface area = mm OL 

Angle subtended by the ith flow tube = os 

Surface area of the ith flow tube at the 

probe surface = garb 


Current flow per flow tube = I/2b = constant 


Current flux in the ith flow tube = T/2b9-r Lb 


Recall that I was previously defined as 





I = “ie , or for this situation, the local current density is 
J 
= 
1 2bo.r Lb 
The local power parameter is G. = ea or after substitution 
Ee 
row 1.2, 
‘ 
ce nis aa G 
i bo. Om 


‘ 
From Figure 16, ¢ = 14°, b = 8 which yields G. = 2.58 G.. 
Applying this correction to the nondimensional temperature-power term 
(0/6) Shifts the theoretical curve so that experimental data more 
closely approximate the theory (Figure 17). 
Since the degree of current density is influenced by probe depth, 
the correction factor must be determined for each experiment. 
Examination of Figure 18 demonstrates the strong end effects on the 


finite cylindrical probe. The cylindrical theory (Equation 11) was 
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Figure 17 
Comparison of experimental data with 
original theoretical curve and 
corrected theoretical curve 
Cylindrical probe L/D = 20 


Data set K 






Corrected theory 


10.0 O 


Original theory 


© Experimental data 
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Figure 18. Liquid crystals demonstrating the 
pronounced edge effects caused by a cylindrical 
probe (L/D = 20) 


48 





based on an infinite cylindrical probe with no end effects. The figures 
show the strong thermal fields at the probe ends developing before the 
effects at the probe center are observed. With time, the conduction of 
heat in the agar tended to reduce the thermal gradient between the probe 
ends and its midpoint, with the resultant growth as shown in Figure 1]. 
A possible explanation for the early observation of thermal fields at 
the probe ends is that radio frequency current density is higher at 
sharp corners than at plane surfaces. The higher current density results 
in greater heating. 

By using a probe with a small surface area, resulting in high 
current density, it was possible. to observe boiling in the agar. Since 
the liquid crystals used in the experiments were active in temperature 
ranges below the boiling temperature of the agar, no datawere taken on 
boiling. However, Figure 19 is a photograph taken of the boiling 
phenomenon. It was of interest to observe formation of a bubble as the 
agar was heated to boiling. When rising, the bubble traveled until it 
encountered a cool region; the bubble then collapsed. This boiling is 
analogous to that occurring in actual surgery. The reason for controlling 
the probe tip temperature is readily apparent when observing the eee 
oration of the jelled agar as bubbles rise randomly towards the surface; 


the analogy in surgery being unpredictable physiological damage. 


49 








Figure 19. Boiling of the agar-water mixture 
caused by excessive probe tip temperature 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The use of liquid crystals for mapping regions surrounding heated 
surgical probes has been clearly demonstrated. When accuracy on the 
order of 1.0°C is sufficient, liquid crystals provide a continuous 
picture of thermal changes as displayed in Figures 20, 21, 22, and 23. 
Through the use of liquid crystals, the spherical theory developed by 
Gengler [5] has been verified. On the other hand, Gengler's cylindrical 
theory has been shown to require further work in order to account for the 
finite cylindrical probe with its strong end effects. 

It is most probable that an actual radio frequency surgical probe 
would be constructed with a length to diameter ratio in the range of 
two to five [1, 2, 15, 16}. Clearly, a theoretical treatment for an 
infinite cylinder would break down at these small length to diameter 
ratios. As a demonstration of this fact, an equivalent surface area 
for a spherical probe was determined based on the surface area of the 
cylindrical probe with a length to diameter ratio of four. Figure 24 
is a comparison of the cylindrical theory, the spherical theory, and 
the experimental data. Relationships between the cylindrical and the 


spherical parameters have been derived based on the fact that the cylin- 


° e & o. 
drical surface area equals the spherical surface area (2mrevql aie Wone) 
Di Pipa D 
(Gn)er; d c fi (Oa ahane J "o_ sphere c ; here 
0 0 


2 | 








Figure 20. Color photograph of liquid crystals in thermal 


field of cylindrical probe (L/D = 4). 
amps /cm? 


Current flux 


0.414 





Figure 21. Color photograph of liquid crystals in thermal 


field of cylindrical probe (L/D = 10). 
amps/cm? 


a2 


Current flux 


0.414 








Figure 22. Color photograph of liquid crystals in thermal 
field of cylindrical probe (L/D = 20). Current flux = 0.414 
amps /cm? 





: 


Figure 23. Color photograph of liquid crystals in thermal 
field of spherical probe 
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Figure 24 
Comparison of cylindrical data (L/D = 4) and cylindrical 
theory with spherical theory and scaled cylindrical data 
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As can be seen, the experimental data is more closely approximated by 
the spherical theory. 

The cylindrical probe with a length to diameter ratio of 20 most 
closely approximated the infinite probe of the cylindrical theory. 
Use of the correction factor to account for current path changes near 
the air-agar interface, tends to bring theory and experiment into 


better agreement. 


B. RECOMMENDATIONS 

Further work, both experimental and theoretical, is needed for the 
cylindrical geometry. Of particular interest are the cylindrical probes 
with length to diameter ratios of five or less since actual probes will 
be in this range. Experimental tnvestigations into the effect of the 
air-agar interface with varying probe depths would demonstrate the 
analogous effects of lesion creation close to the skull; both the air- 
agar interface and the skull are essentially insulators. 

Liquid crystals placed perpendicular to the axis of the cylindrical 
probe and running from the air-agar surface to the bottom of the agar 
could effectively map the distorted thermal region due to the insulated 


air-agar surface. 
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Studies of various tissues should be conducted to fully catalogue 
tissue electrical properties. Such a catalogue would allow the surgeon 
to theoretically determine current strength and duration for a particu- 
lar lesion size prior to the actual surgery. 

Some refinement in the spherical experimentation would reduce 
uncertainty. The major uncertainties were in the thermophysical prop- 
erties of the agar. 

A method of retaining the liquid crystals in the plane of a diameter 
of the probe is required. This is particularly difficult when using the 
cylindrical probe. Placing the substrate in tension in its holder only 
succeeds temporarily. The substrate of mylar eventually stretched and 
sagged. Without having the liquid crystals in the diametral plane, true 
distances cannot be measured for scaling with respect to the probe radius. 

A better potentiometer is reguired for controlling the RFG-2AV gen- 
erator output. The control provided by the manufacturer had excessive 
play. This was partially overcome by careful selection of meter scales 
during experimentation. 

Use of a 35mm camera with a through the lense light meter system is 
suggested for taking data. Most data for this study was recorded using 
the industrial Polaroid film. When used under artificial light, the 
film is very slow. (ASA 24); colors are not well defined which is criti- 
cal when using liquid crystals. The net result with the color Polaroid 
1S a poor reproduction of the liquid crystal colors. Using 35mm black 
and white or tungsten corrected color film gives highly satisfactory 
results when used with the proper light meter settings. Cameras using 
Polaroid film are often moved during film insertion and removal. The 


Slightest movement results in an unfocused picture at the short distances 
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normally used in photographing the probe area. With a 35mm camera, 
there is a minimum of movement during the shutter release and film 
winding processes. The actual time for inserting and removing the 
Polaroid film limits the data that can be taken; this is most true at 
smal] times. Using a 35mm camera, with the fast film advance time, 
allows more frequent data recording. 

A study should be conducted to determine the feasibility of mix- 
ing more than three liquid crystals without a loss of color brilliance. 
The required separation of temperature ranges should be determined 
Since the closer the temperature ranges, the more crystals could 
theoretically be mixed. This would yield more data for each experi- 
mental trial. The temperature ranges, however, must be far enough 


apart so that each color band is displayed clearly. 
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APPENDIX A: RESISTIVITY 


BACKGROUND 

Initially, thermal and electrical properties of the 99.7 percent 
water/0.3 percent agar solution were unknown. It was decided to inves- 
tigate the electrical resistivity of the solution since thermal proper- 
ties could be correlated with the solution's water content [10]. A 
test cell, pictured in Figure 25, was designed and manufactured for 
use in the resistivity investigation. To prevent a standing wave in 
the resistivity test cell, the test cell length could not be an integral 
or a quarter wavelength of the 500 kilohertz Signal used in the experi- 
ments. Using the relationship between frequency (f), wavelength (i), 
and the speed of light in air related to the speed of light in water 
(n); 


8 
= oI meters 
A f 


the wavelength was determined and the test cell designed to prevent 
standing waves. 

A plexiglass tube, 13 cm. in length with 2.5 cm. ID, was used as 
the test cell due to the insulating properties of plexiglass at the 
power levels used. Copper plugs were machined to fit tightly into the 
cylinder with the area of the plug ends known. A BLH two-inch, 
constantan-chromel needle thermocouple was introduced through the wal] 
of the plexiglass cylinder midway between the copper ends. To prevent 
leakage of air into the container, a seal of neoprene rubber was provided 
between the test'cell and the thermocouple shaft. During the experimental 


trials, the test cell was filled with hot agar solution. 
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By carefully capping the test cell, no air was trapped inside. 
The capped cylinder was placed in the generator circuit shown in 
Figure 26. 

Voltage was applied across the circuit for one minute when col- 
lecting data. One minute was chosen since no change in voltage drop 
occurred after this time. Determination of resistivity was accom- 
plished by recording voltages across the circuit and the ten ohm 
resistor along with the distance between the surfaces of the copper 
plugs. The voltage drop across the resistor yielded circuit current 


from Kirchoff's Law: 





The resistance of the cylinder was determined by dividing this current 
into the voltage drop across the cylinder. Resistivity is related to 


resistance by the area to length ratio: 


A 
a. 


It was thus possible to calculate resistivity. 

At the same time the voltage readings were recorded, the needle 
thermocouple was read via the millivolt potentiometer. Data were taken 
as the test cell cooled and again as the cell was heated in a controll- 
able oven (Figure 27). Early experimental trials showed a wide variation 
in agar resistivity (Figure 28). To reduce all plots to a common scale, 


two nondimensional quantities were introduced: 


T - Le 
a <a Temperature 
O 
a ae ) Resistivity 
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MILLIVOLT POTENTIOMETER 
ICE BATH 





RADIO FREQUENCY 
GENERATOR 






RESISTIVITY TEST CELL 


Voltmeter Vy measured voltage drop across the ten ohm resistor 
from which circuit current was determined, 


Voltmeter Vo measured circuit voltage which was used in conjunction 
with circuit current to determine resistivity test cell resistance, 


Resistivity Test Cell Circuitry 
Figure 26 
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Figure 27. The resistivity test cell in the oven 
for the heating phase of resistivity data collection 





Figure 28 
Variation in resistivity values 
as a function of temperature 
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where the subscript indicates the initial conditions of the experimental 
trial. In this way, the plots had common points of origin. 
A linear curve was fit to each set of data points. In dimensional 


form, 


where y was the slope of the curve. Nondimensionalizing yielded 


0 = |-«+ mé 
a 
where m = oa the nondimensional slope. The values of m were 


0 


generated by computer program III (Appendix E). Figure 29 is a sample 
of the nondimensional curves plotted by this program. 

As a result of this portion of the investigation, the great change 
in resistivity over the temperature range of interest (15-60°C) was 
demonstrated. It was determined that the programs written by Gengler [5] 
required modification to reflect the changing resistivity. This was 
accomplished by incorporating the slope of the linear curve from the 


resistivity data into computer programs I and II (Appendix E). 
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APPENDIX B: LIQUID CRYSTALS 


BACKGROUND 

Matter can exist in solid, liquid, or gaseous state depending upon 
the temperature and pressure of the surroundings. In the solid state, 
most compounds are composed of an ordered array of molecules in a 
crystal structure. When heated, the crystal begins to deform as the 
molecules absorb energy and leave the "ground state." The temperature 
at which the deformation begins is known as the melting point [17]; 
it is normally a well-defined characteristic of the compound. 

In 1888 Reinitzer noticed that not all compounds changed directly 
from solid state to an isotropic liquid [18]. The intermediate phase, 
a "“mesophase,' exhibited characteristics of both crystalline structure 
and isotropic liquid. Thus, the name of these substances was coined 
as liquid crystals [18]. Wehnann continued Reinitzer's work in liquid 
cryStals. SPSS Fie the crystals into three categories based on 
molecular orientation; the categories are smectic, nematic, and 
cholesteric [17]. Only the cholesteric crystals will be discussed as 
it is the only type utilized in this report. Interested readers may 
find information on smectic and nematic crystals in the literature 
mre, 1 18,7) 207. 

The cholesteric mesophase has the molecules' long axis parallel and 
coplanar. A model for this structure is a deck of cards, each card 
with one corner bent up such that a slight rotation of the entire 
Stack results [18]. The average displacement of adjacent layers is 


Fifteen minutes of arc per layer [21] with each layer on the order of 
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1000 A thick [22]. The rotation is due to the stretched out non- 
symmetric shape of the crystal molecule [23]. 

The helical rotation of the cholesteric liquid crystal gives rise 
to its optical properties. When transmitting polarized light, the 
helical layers cause the light to be rotated to the left with rotations 
of fifty degrees per millimeter. This property is known as being 
optically active. By way of comparison, quartz, one of the most 
optically active materials in nature, has a rotation of twenty degrees 
per millimeter for polarized light. Another property of the cholesteric 
liquid crystals is crcular dichroism which is a function of temperature, 
material, and the angle of the saGhen light. The incident beam of 
unpolarized white light is split into two components. One component 
is rotated clockwise while the other component is rotated counter- 
Clockwise. The result is that one component is reflected while the 
other is transmitted. To absorb the transmitted beam, a black background 
is required behind the liquid crystals. Circular dichroism is the char- 
acteristic that gives cholesteric liquid crystals their iridescent 
colors. This characteristic should not be confused with the birefringence 
or double defraction of white light by the cholesteric liquid crystals 
(Figure 30). Birefringence is often used as a test for the cholesteric 
mesophase [21]. In the material the incident beam is divided into two 
polarized components which are vibrating at right angles to each other 
and are travelling at different speeds. They are, therefore, refracted 
at different angles and emerge as two beams of perpendicularly polarized 
light. The rotary power of the birefringent cholesteric mesophase is 
as much as 1800 degrees per centimeter when light is transmitted along 


the optical axis normal to a thin sheet [23, 24]. The thickness of the 
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Birefringence: In the material, the incident beam of light is divided 
into two polarized components vibrating at right angles to each other 
and traveling at different speeds. They are, therefore, refracted at 


different angles and emerge as two beams of perpendicularly polarized 
light 





Circular dichroism: The incident beam of unpolarized light is split 
into two components, one rotating clockwise and the other rotating 


counterclockwise. The result is that one is reflected while the other 
is transmitted. 


Figure 30 
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cholesteric film does not alter the predominant wavelength of the 
reflected light as it is rotated; however, the angle of incidence 
and the angle of observation both effect the wavelength of the reflected 
fight. 

The most striking changes in the cholesteric mesophase are those due 
to temperature changes. This mesophase is made temperature sensitive 
by the narrow band of optical reflection which is associated with a 
Change in crystal structure. Cholesteric liquid crystals are also 
sensitive to organic compounds, shear, and electromagnetic radiation 
at high frequencies and power levels. Depending on composition, chol- 
esteric liquid crystals span a range from -20 to 250 degrees C [25]. 
Any specific compound may have its range narrowed to as little as one 
degree or increased to six degrees by varying the compound ingredients 
[26]. If heated far beyond its temperature range and allowed to cool, 
the transition colors occur at lower temperatures, similar to a magnet 
being heated above the Curie temperature [21]. If not thermally shocked 
as just jeeerineet the transition colors remain invariant with tempera- 
ture; the material may be repeatedly recycled. Not all cholesteric | 
liquid crystals react with the same color changes; some display all colors 
of the spectrum starting with red to yellow to green to blue, while others 
may change directly from red to blue. 

Cholesteric liquid crystals will deteriorate in several days if not 
insulated from dust and organic materials. National Cash Register 
Company has been successful with encapsulating the cholesteric liquid 
crystals in spheres 5 to 30 microns in diameter [23]. The result of 
the microencapsulation process is a long life with little deterioration 


from contamination. The NCR process reduces the angle dependence of the 
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crystal as far as color is concerned although resolution is degraded 
and response time is increased approximately ten percent [22]. The 
microencapsulated crystals are used in industry since they are not 
readily contaminated; they were used in this experiment. 

Liquid crystals have many engineering applications. In temperature 
measurement, cholesteric liquid crystals are primarily used. Guidelines 
for use of liquid crystals in thermal mapping may be summarized as 
follows: 

a) The heat capacity of the object should be greater than that 

of the liquid crystals; 
b 


ee? 


The size of the object must be great enough to resolve the 

colors of the liquid crystals; 

c) The rate of temperature change must be slow enough for the 
liquid crystal to follow. (A maximum time of about 0.2 seconds.); 

d) The temperature range under study must be within the range of the 
liquid crystal; 

e) The surface of the object must be capable of holding liquid 
crystals unless an adhesive backing is used; 

f) The surface must be black to display the liquid crystals; 

g) The temperature pattern must be correlated with another 

material property such as voids present in nondestructive 


testing or power dissipated per unit area. 
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APPLICATION AND CALIBRATION 

Reference 27 lists many techniques for application of liquid 
crystals; in this investigation, a commercial artist's air brush was 
utilized. Three crystals, designated R-27, R-37, and R-53, were mixed 
equally. After standing overnight, the cholesteric substances tended 
to separate from excess water. The water was carefully removed leaving 
a more concentrated liquid crystal which resulted in a more brilliant 
color display. The concentrated mixture was sprayed onto the transparent 
substrate in three thin coats, allowing adequate drying time for each 
coat. To prevent any absorption of moisture from the environment, dry- 
ing was conducted under a 40 watt lamp. After the three coats of liquid 
crystals were dried, two coats of black paint were applied. Testors 
Flat Black Spray Enamel was used since it contained no chemical which 
might have reacted with the crystals to cause a shift in the color- 
temperature relationship. The black paint was required to absorb the 
light transmitted through the liquid crystal, as only a small portion of 
the light was reflected. To provide a watertight system, a layer of 
polyurathane varnish was applied over the black paint and along each edge 
of the substrate. Thus, the crystals were sealed from all contaminants. 

A second type of liquid crystal, designated ROCHROME, was also used 
during the experiments. Rochrome, manufactured by Hoffman-LaRoche, 
consists of liquid crystals sandwiched between a layer of mylar and a 
black adhesive backing. Rochrome is cut to the desired Shape and applied 
directly to the object under study. Again, the edges were sealed to pre- 
vent water contamination and separation of the Rochrome from the surface. 
The surfaces to which the Rochrome and the liquid crystals were applied 


in this investigation were mylar and glass. The mylar was used with the 
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probe apparatus while the glass was used during the calibration 
procedure. 

Calibration was conducted using the Rosemount Engineering Company 
Variable Temperature Bath (Model 913A) in conjunction with its 
commutating bridge and platinum resistance standard thermometer. 

National Cash Register liquid crystals (R-27, R-37, and R-53) 
were prepared by spraying them onto a glass slide, coating with 
black paint, and sealing with polyurathane varnish. Rochrome was 
prepared by sticking a small piece to the slide and sealing the edges 
with polyurathane. The slides were placed in the bath at approximately 
one-third of the bath's depth and close to the center of the bath 
diameter. Effects of the bath container edges were thus minimized. 
The bath temperature was raised slowly as the slides were observed. 
To minimize variations between observers, the same observer read the 
Slide colors while an assistant operated the equipment during all 
calibration procedures. | 

As the fan aeeetire was raised, the liquid crystals reacted, each 
at its own temperature level. When the slide had turned one color, 
1.e. red to green, the commutating bridge reading was recorded; then 
the temperature of the bath was increased to determine the next color 
change. A conversion table was used to convert the bridge readings to 
temperature. The uncertainty in this temperature-color calibration 
procedure was calculated to be +0.5°C which agrees with Reference 27. 
Table B-1 lists the crystals used and their calibration temperature- 


color relationships. 


ae 





Crystal 
R-27 
Rochrome 
R-37 
R-53 


Red 


2 Gres eee 
Zoo 
Sane 
SOSoMe 


Table B-] 
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Green 
27.9°C 
31.5°C 
37.0°C 
52.1°C 


Blue 


28.9°C 


34.2°C 
g6,5—C 
a3 JSC 





Apeevome, C: PLUK PLOTS 


Gengler's theoretical analysis of the cylindrical probe modeled the 
probe as an infinite cylinder with negligible end effects; the current 
flow should have been a one-dimensional problem. With the experimental 
cylindrical probe, end effects were pronounced, prompting an investiga- 
e7on oF Current paths in the probe vicinity. The area of interest was 
within ten probe diameters. The problem was approached in two phases. 
First, the horizontal plane was examined using TRUMP [28, 29], a computer 
code, while in the second phase, teledeltus paper was used to depict the 
isopotentials in the vertical plane. 

TRUMP solved the general partial differential equation Ve = 0. 

It was postulated that the radio frequency potential steady state 
condition was established instantaneously with power application. 

A nodal network was constructed for use with TRUMP (Figure 31). Due 

to symmetry of the probe-cell system, only one quadrant of the system 
was modeled. Three cylindrical probes, with length to diameter ratios 
of 4, 10, and 20 were used with mo sets of boundary conditions. The 
first boundary conditions approximated a small finite cylinder inside 
an infinite cylinder by the use of two conductive sides of the nodal 
boundaries for the outer cylinder wall and two insulated sides of the 
nodal boundaries for the endcaps of the outer cylinder. Figures 32, 
33, and 34 are plots of the isopotentials and lines of constant current 
for each probe length to diameter ratio. Inspection of these figures 
in the area within ten probe diameters, shows the current lines orthogonal 


to the probe surface. It was, therefore, expected that the current flow 
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and the resultant thermal problem would be one dimensional. For 
comparison, a second set of boundary conditions surrounded the probe 
with conductive surfaces. Figures 35, 36, and 3/ display the iso- 
potentials and constant current lines for this situation. 

The second portion of current path investigation was undertaken 
when experimental data did not agree with the theory. As stated earlier, 
in this second phase, isopotentials in the vertical plane were examined 
using the same sets of boundary conditions as in the first phase. 

The probe-cell system was scaled on teledeltus paper. Boundary 
conditions were simulated using conductive silver paint and metallic 
tape. Figures 38 and 39 are plots of the isopotentials and constant 
current lines under the two sets of boundary: conditions. It is seen 
that the insulated agar surface caused current to bend towards the 
conductive boundary. Since heating is a function of current per unit 
area, the higher current density resulted in higher heating effects 
in the experimental model. The higher heating effects were the major 
cause for the difference between theory and experiment with the 


cylindrical probe. 
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SH Teledeltus plot of 
OX isopotentials and constant 


\S current lines in a plane 
perpendicular to the cyl- 
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Figure 39 
Teledeltus plot of 
isopotentials and con- 
stant current lines in 
a plane perpendicular 
to the cylindrical 
probe axis. 


Free surface of agar- 
water mixture acting 
as conductive boundary. 
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APPENDIX D: SAMPLE UNCERTAINTY ANALYSIS 


Reference 30 categorizes error in experimental investigation as 
fixed error, variable error, and mistakes. Fixed error is defined 
as a constant error such as equipment consistently reading low. 
Variable error is caused by variations in sensors, an example of 
which is having two different observers taking readings in separate 
experimental trials. Finally, mistakes may occur during experiments 
by reading the wrong scale, etc. If it is assumed that the error, 
j.e. uncertainty, of each primary variable has the same probability 
of occurrence, it is then possible to calculate the overall uncer- 
tainty in each value. As an example, define P, the result of some 
experiment as, 

pP 7= P(x, 5 Xos Xa5+++X)) 
If w. is the uncertainty in P, ane W1s Wos Wzs ++-W, are the 


p n 
uncertainty in each primary independent variable, then the relative 


uncertainty in P is defined as 





In the experimental trials on the spherical probe, the quantities 


calculated were 


125 V25 
meee lene 2 
O T ro 2 1 ro k of 10 
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Re ss =&4) ome 


Deriving the expressions for the relative uncertainty in each 


quantity 
1/2 
eG He Ne oTe\2 Zu \2 
6... mye O + 
G - K T R 
O O € 
no 2 
8 
ae 
R 
1/2 





For all independent variables except resistivity 0. it was possible 
to estimate the uncertainty either directly from observing equipment 
(voltage, current) or from available tables (ao, k). For resistivity 


it was necessary to calculate uncertainty using 


V a 


OO SS HE) 


Following the example of Equation D-1, 


85 








G) Uy) 2 WW 
EO 3 Vect i f veer 
ee! li V V 
Po Corea ik cee 
Z Z Z 
WR an 7 V2 
10 D L 
R | ie 
10 


Table D-] lists values of uncertainty for each independent variable 
and the value of the variable used in the calculations. The following 


values were calculated for uncertainty: 


0 
—2° = 0.036 (3.6%) 
Po 
ae 
= = O.1217 (12.22) 
0 
“9 
x = 0.019 (1.9%) 
YR 
= = 0.029 (2.9%) 
G) 


La 


= OmlOZ (10e27) 


= 


To obtain uncertainty in the quantity 6/G the following relationship 


was used 
”9/G w\2 (°G \2) 1/2 


Bee ee ee 
6/G 0 G 


0 


Using the values above, 


“0/G_ 
“0/6 =e Omlicor lice. 
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As may be seen from the last column of Table D-1, the largest 
uncertainties arise from the thermophysical properties, i.e. thermal 


conductivity (k), thermal diffusivity (a), and resistivity (°o). 


Table D-1 

Quantity i Uncertainty “i pewesal Uncertainty 
eet = 80.0 volts 0.0025 volts 0 
Vp = 0.118 volts 0.0025 volts 0.02 
D = 2.540 cm 0.025 cm 0 
L = 11.95 cm 0.025 cm 0 
Rig = 10.0 ohms 0.3 ohms 0.03 
V = 0.60 volts 0.0025 volts 0 
-— = 0.241 cm 0.00] cm 0 
i. = 2320= C OnZor aC 0.01 
it = 37.0 °C OnNGeee 0 
- * 0.380 cm 0.01 cm Os OZ 
ro = 0.90 cm 0.0] cm 0.0] 
t = 30.0 sec 0.025 sec 0 
k = 0.628 W/m°K 0.0628 W/m°K Oo 
O = 0.001512 m¢/sec 0.0001512 m¢/sec 10.0 
p = 3250.0 ohm-cm 117.0 ohm-cm onbU 
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